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1. Motivation 

Theories different from QCD have seen a considerable amount of recent attention from the 
lattice gauge theory community. In particular, the properties of SU {Nc) Yang -Mills gauge theories 
with Nf light fermions have been the subject of numerous recent lattice studies [1, 2, 3, 4, 5] in 
the context of searching for novel non-perturbative physics which could play a role in building 
strongly-coupled models of physics beyond the standard model, e.g. technicolor. We will consider 
here only the case Nc = 3, with Nf fermions in the fundamental representation of the gauge group. 

The infrared dynamics within this class of theories is known to undergo a transition as Nf is 
varied within the range <Nf < N^ for which asymptotic freedom is preserved. It is well known 
that for Nf just below N"\ the jS -function of the theory shows a weak infrared (IR) fixed point 
[6, 7], leading to conformal behavior in the IR. This is markedly different from the confining IR 
behavior of a QCD-like theory (small Nf); we expect a phase transition at some point as Nf is 
varied between the two extremes. Thus far, lattice studies have supported the picture that there is 
a "conformal window" Nj < Nf < N"^, with the fixed-point coupling strength becoming strong as 
Nf approaches the critical transition value Nj. 

Lattice simulations are necessarily performed at finite fermion mass nif, and must then be ex- 
trapolated to the physical point or the chiral limit. Thus, chiral perturbation theory ix^T) is an in- 
valuable tool in connecting our simulation results at various Nf with continuum physics. Although 
most interest thus far has been on SU (2) and SU (3) ;^PT for application to QCD, generalization 
of ;tPT to the breaking of an SU (Nf) x SU (Nf) chiral symmetry is quite straightforward, with the 
Nf counting factors well-known up to next-to-leading order (NLO). However, the values of the 
low-energy constants appearing in the chiral Lagrangian are not identical between the SU{2) and 
SU{Nf). 

Of particular interest is the magnitude of the fermion condensate (yy) relative to the Gold- 
stone decay constant F, encapsulated in the ratio (y7i//)/f ^. In extended technicolor models, there 
is significant tension between the requirements of evading experimental bounds on flavor-changing 
neutral currents (requiring a large ultraviolet cutoff Aetc)> and generating standard-model quark 
masses which go as lJifY)Tc/ ^etc^ where {yv)tc is the technifermion condensate cut off at the 
scale Aetc- A phenomenon known as "walking" can significantly enhance the ratio (yw) /^^ rel- 
ative to QCD [8, 9, 10], allowing generation of the correct masses without violating experimental 
bounds (see e.g. [1 1] for a brief review.) Continuum studies based on Feynman graphs have sug- 
gested just such a significant enhancement for theories with many flavors Nf, mar (but just below) 
the transition value A^^ [12, 13]. 

Determining the evolution of (yv) /F^ is a primary goal of the Lattice Strong Dynamics 
(LSD) collaboration. Initial results were described in Ref. [14], comparing the well-known Nf = 2 
theory to the case of Nf = 6, concluding that there is significant enhancement at Nf = 6. In this 
supplement to Ref. [14], to appear in the Proceedings of 6th International Workshop on Chiral 
Dynamics, we describe more completely the use of ;^PT and the extrapolation to m = 0, with 
particular focus on the difficulties that arise at Nf = 6. 

In section 2, we review the ;^PT formulas for general Nf. In section 3, we detail our simulation 
methods and parameters. In section 4, we perform a variety of chiral fits to our data at both Nf = 2 
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and Nf = 6, with comparison between the two cases. In section 5, we make some concluding 
remarks. 



2. ;t:PT at general A^/ 



The appUcation of ;^PT to theories with a general number of light fermion flavors Nf (all 
with mass m) is straightforward. Aside from the unknown Nf dependence of the low-energy con- 
stants themselves, Nf appears as a loop counting factor in NLO formulas and beyond. The general 
formulas at NLO for A/^, and {YV)m are [15]: 
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where a,- = 8(4;r)2L; and t], = 8(4;r)2//,- are just the conventional low- and high-energy constants 
of X?T, rescaled to values of 0(1). Both the analytic terms and the chiral logarithms scale with the 
number of fermions Nf; except for the log correction to , the size of the NLO terms increases 
relative to the leading order as Nf increases. The lone high-energy constant r\2 which appears in 
{\j7\l/)m includes a "contact term", a quadratic divergence in the ultraviolet cutoff (here I /a, where 
a is the lattice spacing). We therefore expect the analytic term linear in m to dominate the chiral 
expansion of {\j7Y)m- 

Without fitting to additional observables or resorting to partial quenching, we cannot distin- 
guish the individual low- and high-energy constants in eqs. (2.1) to (2.3) above. Therefore in our 
analysis, we adopt a simpler notation by combining many of the coefficients: 
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The quantity zm is the expansion parameter of ;^PT, so that 1 /z gives a rough estimate of the mass 
scale at which we expect perturbation theory to break down. 

At next-to-next-to-leading order (NNLO), six additional coefficients appeai^, corresponding 
to the and log m terms in the chiral expansion. The coefficients of the leading non-analytic 
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terms m^(logm)^, like the mlogm terms above, are completely determined by Nf, and have recently 
been computed for general Nf [16]. Extending the notation of eqs. (2.1) to (2.3), we write 

l2 
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Because of the dominance of the linear term in the expansion of {\j7\i/),„, when fitting our data to 
eqs. (2.8) to (2.10) we will fix the unimportant NNLO analytic terms ac20 = otcii = 0, in order to 
obtain a fit with more degrees of freedom. 



3. Simulation details 

The study of unexplored aspects of chiral dynamics is a difficult proposition on the lattice, 
as lattice fermion discretization typically breaks chiral symmetry explicitly. For this study, we 
employ domain wall fermions [17, 18], which have good flavor symmetry properties and break 
chiral symmetry only by an exponentially small factor. In addition to the domain wall fermion 
action, we use the Iwasaki improved gauge action [19]. 

Gauge configurations are generated using the CPS application package, part of the US QCD 
software library. Evolution in configuration space is performed via the hybrid Monte Carlo (HMC) 
method, optimized with a three-level symplectic integrator, a single level of Hasenbusch precondi- 
tioning, and chronological inversion. 

The lattice volume is set to 32^ x 64. For the length of the fifth dimension we take = 16, 
and the domain-wall height is set to mo = 1.8. Input fermion masses are varied from mf = 0.005 to 
0.03. Raw data are blocked over sets of 50 trajectories before analysis, in order to reduce the effects 
of autocorrelations; our runs are not long enough to perform a complete analysis of autocorrelation 
times. 

Since we do not extrapolate Ls — )• °°, there is some residual chiral symmetry violation, which 
can be encapsulated in a residual mass m^es- The magnitude of the residual mass in the chiral limit 
is rures = 2.60 X 10"^ in the Nf = 2 theory, and m^es = 8.23 x 10""^ for Nf = 6. To take the chiral 
limit, we extrapolate in the total fermion mass m = mf + m^es- The variation of m^^ over the range 
of m/ used is extremely small compared to m, and we neglect it here. 

In order to obtain a correct sampling of the partition function in finite volume, evolution of 
global topological charge Q is important. Our lightest evolutions at mj = 0.005 show signs of very 
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slow evolution in topological charge, which can lead to significant systematic shifts, particularly 
in {\j7Y)m and [20]. Due to this unknown systematic effect, we exclude the mf = 0.005 results 
wherever possible below, and caution that any fits which do contain data at this mass point ai^e 
questionable. 

4. Chiral extrapolation 

We will consider an assortment of chiral fits, at both NLO and NNLO, varying both Nf and 
the mass ranges included in the fit. The set of fits to be considered is detailed in Table 1 . Resulting 
best-fit parameters are shown in Table 2. 

We begin at Nf = 2 with fit type A, which is simply an NLO fit to all of the available data 
(excluding mf = 0.005, for the reasons given above.) The ;^^/d.o.f. is 6.5 with 4 degrees of free- 
dom; given that our error bai^s are most likely underestimated due to the relatively small number 
of trajectories gathered, this indicates a reasonably good fit to the data. From Table 2 we note that 
aM and are 0(1) numbers, while Oc is very large so that the linear term dominates the {\j/W)m 
extrapolation, as expected from the presence of the "contact term". We expect the breakdown of 
chiral perturbation theory for zm f« 1; since 1/z = 0.036, the expansion should be good over the fit 
range used. 

Despite our concerns about the nif = 0.005 data, a natural question at this point is whether 
the NLO fit can be extended to include them. We attempt fit type B in order to investigate. The 
X^/d.o.f. rises significantly to 36.2, now with 7 degrees of freedom. The inability to incorporate 
these lighter-mass points into the original chiral fit is most likely due to the presence of large 
systematic errors in the data which are not reflected in the error bars, as we anticipated. 

Attempting the same basic NLO fit (type A) a.tNf = 6 is the logical next step, but unfortunately, 
the data does not support such a fit. As shown in Table 2, the values of the best-fit parameters are 
fairly close to those from the Nf = 2 fit, but the ;^^/d.o.f. is now 50.5, indicating a very poor fit 
to the data. Both fits of type A are plotted in Figure 1, and the poor quality of the A^j = 6 fit 
is immediately obvious. There is clear tension in the fit between the intercept values F and B 
expected from more naive extrapolation of the data, and the magnitude of the NLO terms governed 
by z B/F^. Including the trif = 0.005 data (type B) does not improve the situation, with roughly 
the same best-fit parameters and a somewhat larger ;|^^/d.o.f.. 



Fit label 


Order 


mf range (M^ ) 


mf range (Fm) 


mf range i{\j7¥)m) 


Nclof 


A 


NLO 


0.01-0.02 


0.01-0.02 


0.01-0.02 


4 


B 


NLO 


0.005-0.02 


0.005-0.02 


0.005-0.02 


7 


C 


NLO 


0.01-0.02 


0.01 


0.01 





D 


NNLO 


0.005-0.02 


0.005-0.02 


0.005-0.02 


3 


E 


NNLO 


0.01-0.03 


0.01-0.03 


0.01-0.03 


6 



Table 1: Types of chiral fit considered. Both the order of the fit and mass ranges of data included for M,^,, 
Fi„ and {\j/\ir),„ data are varied. Fits A through D are performed for both N/ = 2 and Nf = 6; fit E can be 
performed only at Nf = 6, where data at = 0.025, 0.030 are available. 
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m m 

Figure 1: Chiral extrapolation of the quantities M^/2m and F„, with ;|fPT fits of type A. Data points used 
in deriving the fits shown are dark red squares (Nf = 2) and dark blue triangles {Nf = 6); data excluded from 
the fit are lightly shaded. Best-fit results with l-a error bands are shown in corresponding colors. 




Figure 2: Chiral extrapolation of the quantities M^j2m and f],,, with jPT fits of type C at both Nf = 2 
(red) and 6 (blue). Data symbols and colors are as in Figure 1. 



The failure of NLO ;^PT at Nf = 6, over the same range of which yielded a good fit at 
Nf = 2, is not entirely surprising in retrospect. Most of the NLO terms in eqs. (2.1) to (2.3) scale 
linearly with Nf, so that if we keep z fixed, then as we go from Nf = 2 to 6, we expect to have to 
work in a mass range which is a factor of 3 smaller to keep the size of the NLO terms fixed relative 
toLO. 

An exception is the slope of the Goldstone mass squared, eq. (2.1). There, the chiral log is 
actually suppressed by l/Nf (although there is still a linear increase with Nf buried in the linear 
coefficient Um-) Indeed, by inspection the variation of M^„/2m over the mass range O.Ol < nif < 
0.02 is quite small, both at Nf = 2 and Nf = 6. We might therefore expect that a chiral fit to just 
the Goldstone mass data might be more successful. This is fit type C, which includes the range of 
Ml data 0.01 < Mf < 0.02, along with the lightest data point nij = 0.01 for both F,„ and {\I/V)m 
in order to fix ap and ac- We have zero degrees of freedom, so evaluation of goodness of fit is 
difficuh. 

Results for fits of type C are shown in Figure 2. Both fits of type C lead to a large value of 
the parameter z, pointing to an early breakdown of the chiral expansion; in both cases, the NLO 
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Figure 3: NNLO chiral extrapolation of the quantities M^j2m and F„j, with fits of type D at Nf = 2 and 
type E at Nf ~ 6. Data symbols and colors are as in Figure 1. The ratios of NLO (red: Nf = 2, blue: Nf ~ 6) 
and NNLO (orange: Nf — 2, green: Nf = 6) terms to LO in the F„, extrapolation are also displayed, showing 
that all terms become quite large relative to LO near the bottom of the mass range used. 

terms in f„, are clearly large compared to the leading-order value at quite small m, with fine-tuned 
cancellation of the NLO terms occurring in order to match the single point contained in the fits. 
Although we have no ;^^/d.o.f., the type C fit curves clearly do not match the Fm data except at the 
point nif = 0.010 included in the fit. 

Finally, we attempt the NNLO fits of eqs. (2.8) to (2.10). As discussed we fix ac20 = cxcii = 
0, so that there are 9 free parameters. In order to obtain a constrained fit, we must include the 
mf = 0.005 data at Nf = 2 (fit type D.) At Nf = 6, we can use the range 0.01 < m/ < 0.03 (fit type 
E), which we expect to yield a better result due to the systematic effects in the WJ/ = 0.005 data. 
Results are plotted in Figure 3. Fit type D, much like type B, includes the spurious nif = 0.005 
data, so the resulting large values of ;^^/d.o.f. are not surprising. Fit type E at Nf = 6 yields an 
apparently good fit to the data based on ;^^/d.o.f.. However, separating out the NLO and NNLO 
terms relative to the leading order shows that both terms are equal to or larger than the leading 
order even at the bottom of the mass range being fit to (Figure 3). In addition, the analytic terms 
are generally large and strongly covariant. The apparent goodness of fit is thus inconsistent with a 
convergent chiral expansion, despite the value of ;^^/d.o.f.. 

5. Conclusion 

Based on initial results from the LSD collaboration, we conclude that NLO ;^PT gives a self- 
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consistent fit to our Nf = 2 data, but does not yield a satisfactory fit to the data at Nf = 6 for our 
current mass range 0.01 <mf< 0.02. Alternative fit ranges and the inclusion of NNLO terms also 
fail to give an acceptable Nf = 6 fit. As noted above, due to the linear scaling with Nf of several 
terms we might expect that at lighter mass points (by a factor of 3) will be necessary in going from 
Nf = 2toNf = 6. The inclusion of known finite- volume and finite-topology terms [20] in the chiral 
fits should allow us to fit to lighter masses without systematic errors and without greatly increasing 
the computational cost. Also, simulation at additional masses, particularly Ughter ones, may allow 
for a better NNLO chiral fit which is consistent with a good chiral expansion. 

Alternative approaches, including partially quenched analysis or the use of eigenvalue meth- 
ods, may provide another way forward for the extraction of chiral quantities at large Nf, and will 
be explored by our collaboration in future studies. Detailed analysis of the chiral extrapolation for 
the quantity {\j/\j/)/F^ and of the physical implications of our results are carried out in [14]. 

We thank LLNL and the Multiprogrammatic and Institutional Computing program for time 
on the BlueGene/L supercomputer. This work was supported partially by DOE grants DE-FG02- 
92ER-40704 (T.A., E.N.), DE-FG02-91ER40676 and DE-FC02-06ER41440, and NSF grants OCI- 
0749300, DGE-0221680, PHY-0427646 (A.A., R.B., R.C.B., M.A.C., S.D.C., C.R., D.S.), PHY- 
0835713 (M.A.C), and PHY-0801068 (G.E). It was also supported by the DOE through ANL 
under contract DE-AC02-06CH11357 (J.O.), and the DOE Office of High Energy Physics through 
LLNL under contract DE-AC52-07NA27344 (M.C., PV.). 



Nf 


label 


7 


F 






ccc 


zVd.o.f. 


2 


A 


28(16) 


0.0209(4) 


0.31(62) 


0.64(47) 


83(29) 


6.50 


6 


A 


25(11) 


0.0188(36) 


2.5(1.4) 


0.1(1.1) 


194(24) 


50.5 


2 


B 


44(9) 


0.0184(13) 


-0.10(11) 


0.90(15) 


58(6) 


36.2 


6 


B 


27(10) 


0.0179(29) 


2.4(1.2) 


0.20(96) 


204(21) 


74.2 


2 


C 


77(12) 


0.0171(5) 


-0.50(10) 


0.74(11) 


25(10) 


OO 


6 


C 


138(63) 


0.0133(14) 


-0.26(10) 


2.1(1.1) 


28(10) 


oo 


2 


D 


16(8) 


0.0259(39) 


20(18) 


-6.6(5.9) 


119(45) 


30.0 


6 


D 


17(13) 


0.0214(66) 


25(37) 


-14(14) 


237(73) 


14.9 


6 


E 


17.9(5.5) 


0.0217(28) 


1.8(7.6) 


-2.2(4.0) 


206(20) 


6.08 


Nf 


label 


«M20 


CCM21 


CCF20 


0CF21 


ac2o 


ac2i 


2 


D 


32(77) 


81(104) 


-5(17) 


-24(29) 


O(-) 


O(-) 


6 


D 


21(83) 


92(167) 


-7(41) 


-51(74) 


O(-) 


O(-) 


6 


E 


4.6(1.3) 


14(18) 


1.7(0.8) 


-10(11) 


O(-) 


O(-) 



Table 2: Chirally extrapolated quantities and fit parameters, based on the assorted ;i;PT fits considered. 
All fits shown are joint fits between the three quantities {\j/\i/)m, F,„ and M,^,. NNLO analytic coefficients 
ac20, Cicii are fixed to zero, due to the dominance of Oc- 
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